A note on the quantum-tail effect on fusion reaction rate 
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A study is made of the power-law tail effect in the quantum particle distribution 
over momentum on the nuclear fusion reactions. Our results do not support the 
idea of averaging the fusion reaction cross-section over the momentum distribu- 
tion postulated and used in many publications. 
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While in classical statistics all systems of particles in equilibrium have single-particle mo- 
mentum distribution (MD), n(p), of Maxwell-Boltzmann form, the MD in quantum statistics 
has a non-Maxwellian form, contains a non-exponential tail [1-4] and plays a role central to 
our understanding of interacting quantum particles systems. n{p) has to satisfy the sum 
rules 

J n(p)d 3 p = 1, (1) 
^ J n{p)p 2 d 3 p = T k , (2) 

where Tk is the kinetic energy. 

In order for the kinetic energy integral, Eq.(2), to remain finite, the MD should decline 
faster than p~ 5 . The large p-behavior of n(p) has been considered in many papers [1-4]. 
It was shown that the interaction between quantum particles leads to the appearance of a 
power-law tail in the MD. It was found for the first time in Ref.[2] that the zero temperature 
MD for an interacting electron gas should go as 1/p 8 for large p. 

A general expression for n{p) can be written as [5] 

n(p)= J f 7 (E,p)dE, (3) 

where the energy-momentum distribution f^(E,p) can be represented as [5] 

f 1 (E,p) = f(E)5 1 (E,p), (4) 

f(E) is the occupation number (Fermi-Dirac, Bose-Einstein or Maxwell-Boltzmann), 8 y (E,p) 
is the spectral function [5] 

l(E,p) 

6 " {E ^ = *[(E - e„ - WW + l 2 (E,p)} ' (5) 

e p = p 2 /(2m), and ^y(E,p) and A(E,p) are solutions of "complex set of integral equations" 
[5]. 

The finite value of ^{E,p) leads to the appearance of power-law tails in the MD. 
After integration over momentum of f y (E,p), it is easy to obtain that the energy distri- 
bution 

n E {E) = Jf 7 (E,p)d 3 p (6) 

remains exponential [6]. The difference between the two distributions, Eq.(3) and Eq.(6) is 
related to the quantum uncertainty. 



Refs.[6,7] have been suggested to average reaction rates over n{p) rather than distribution 
over energy n E (E) (see also Ref.[8]). Since (i) in this case quantum tails might produce dra- 
matic effects on the rates of nuclear and other reactions in a medium and (ii) the provocative 
suggestion [6,7] has been used in a many papers [6,7,9-18], it is the purpose of present work 
to examine the validity of the approximation [6,7]. 

We stress here that directly measured cross sections of low-energy nuclear reactions in a 
medium, when available, are higher than the expected values [19-32]. To date, the observed 
enhancement factors are not completely understood. An example of the effort where envi- 
ronment considerations have been carried out for low-energy processes are the Refs. [33-35] 
(see also Refs. [36,37] where the tunneling of bound two-body systems through a potential 
barrier have been considered). 

Let us first study the effect of the binding of the deuteron inside the hydrogen atom on 
the the example of the p — d fusion reaction. The bound deuteron does not have a definite 
velocity. From its wave function one can calculate the deuteron velocity distribution 

8 1 

n(vd) = -0-3 /-, ; / — / x 2 x 4 > (7) 
n 2 vl (1 + (v d /v ) 2 ) 4 

where vq = (m e /md)e 2 /K, and m e and rrid are the electron mass and deuteron mas, respec- 
tively. One has 

J n(v d )d 3 v d = 1. (8) 
Now, the bound deuteron and the proton have a relative velocity 

Vrel = \Vd - Up I (9) 

and < av > is then 



f 1 E 

<av>= J S{E)—exp(-7iJ -g-)v re in{v d )d 3 v d) (10) 

where E = /iv 2 el /2, \i is the reduced proton deuteron mass, S(E) is the astrophysical ^-factor 
and E G = 2e 4 fi/ti 2 . 

In the case of small v p (v p < v ) the < av > value, Eq.(10) corresponds to the screening 
energy E scr xs 300eV. This result is not in agreement with three body adiabatic calculations 
[8,38] (more than 10 times larger). In our future work we will consider an application of the 
Sturm-function method in the formalism of the three-body Faddeev-Hahn equation [38,39] 
for this problem. 
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We consider a system of N identical quantum particles carrying the unit positive charge, 
e, and contained in a volume Q (periodic boundary conditions) with an uniform external 
background field of the opposite sign which neutralizes the total charge of the system. The 
time independent n-th state wave functions of the system $ n (fi, r 2 , ...f^) with the energy 
E n are assumed to be normalized. 

At a thermal equilibrium, n{p) can be written as 



T,nfn(P, 



-E n /(k B T) 



<P) = ^^-E^T) » (H) 



where f n (p) is the probability to find a particle with the momentum p in the n-th state <3> n 

N 

- I pM r_r )<b*< r tv, -r^d) (r 1 tv. r^Arl 3 rrl 3 r' 

(27T) 



If 

/n(p) = 7^ / et^)$;(f,f 2 ,...r1v)$n(^y2,...nv)^Vn^. (12) 

l 27r J J i=2 

A generalization of the Kimball method, Ref.[l], leads to the following large-p behavior of 

Up) 

9 N — 1 rr? 2 p 4 

Bm/ B (^) = W 0)| 2 -— ft^, (13) 

where 

AT 

|^ n (0)| 2 = fi/ |$„(f,f,f 3 ,...r-W)|Vrn^ (14) 

J i=3 

and m is the particle mass. Substituting the large p asymptotic, Eq.(13), into Eq.(ll) we 
find the large momentum tale of the momentum distribution n{p) at the temperature T 

lim n(p) = -hp— ^\V(0)\ 2 , (15) 

where 

is the contact probability of finding two particles at zero separation, and p = (N — 1)/Q is 
the density for N » 1. 

We thus see that the large-p behavior is governed by the contact probability of finding 
two particles at short distance. Although, the power-law tail in the momentum distribution 
was observed in many papers (see, for example [1-4]), it was shown for the first time in Ref.[l] 
that for the ground state coefficient of the 1/p 8 term is proportional to the zero separation 
probability |-?/>o(0)| 2 - 
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Now, we consider nuclear reaction between nuclei, i and j, under conditions that exist in 
stellar interiors. The nuclear cross-section, for small collision speed is approximated by 

a^E) = ^lexp(-n^jE) « ^ exp(-7r^7s), (17) 

where a slowly varying function with E (E << Eg) Sij(E) is called the astrophysical S- 
factor, Eq = 2p i jZfZ'je' l /fi 2 is the Gamov energy, //^ = mirrij/(mi + rrij) and Z denote 
charge numbers. The last equation (17) is practically exact - for example, the error in the 
replacement S(E) by 5(0) for the proton-proton fusion, p + p -^ 2 D + e + + v e is only about 
0.5 % [41] for energies corresponding to temperatures in the center of the Sun. 

The number of reactions between nuclei of i and j species at number densities pi and pj 
with a relative kinetic energy E is calculated as 

= j^h a - u ' :)r ^ (18) 

where Vij = ^2E / pij. 

The Gamov rate is calculated by averaging Eq.(18) with the Maxwell-Boltzmann distri- 
bution fMs{E) at temperature T « Ec/hs 



'"-^-^te^-isf*- (19) 

The result yields the Gamov rate [42] 

Rfj = < >m-b= J f M ~B(E)cr l3 (E)v l3 dE. (20) 

It is possible to rewrite the rate in terms of the contact probability (the square of the wave 
function at the origin [43]) 

R^E) = hS iS (E) P \ P f 77p2 hM0)| 2 , (21) 

where ip^ is the Coulomb wave function with a normalization such that 

' |^-(r)|Vr = ii (22) 



/ 

Jn 



in 

over a large volume fl [43] . Indeed, the square of the wave function at the origin then takes 
on value 



2 KyjE G /E 



hMo)l 2 = - h= '=- — • (23) 

exp(nJE G /E) - 1 



For the iV-body system, Eqs. (11-16), the number of binary fusion reactions per unit time 
and unit volume in the n-th state § n (ri,r 2 , •••'Hv) is 

Rn = S{m— 2 \i>nm 2 = 5(0)^^3 Jim (p 8 f n (p))- (24) 

In (24) |*0 n (O)| 2 and f n (p) are given by (14) and (13), respectively. 
At the thermal equilibrium 

Y^ n e -E n /(k B T) v ^2(me 2 ) 3 v 11 v 1 

Now, we want to see if the rate R, Eq.(25), and the rate Rq, calculated in Ref. [13] by 
averaging the pp-fusion reaction cross-section over quantum momentum distribution [3], 
stand in contradiction to each other for the pp fusion in the star. To do so we note that, in 
the Galitskii-Yakimets approximation for n{p) at large momentum [3], ratio R/Rq is 

R n 

* " 10 ' 

where the standard solar model parameters [44] are used. 

Clearly, the error in the averaging the fusion reaction cross-section over quantum mo- 
mentum distribution is that one has neglected the coupling between the various probability 
amplitudes of velocity which is introduced by the quantum uncertainty. 

In conclusion, we summarize the main points of this letter. 

(i) We have considered nuclear motion inside the atom and have found that for the p — d 
fusion in the case of small v p , the < av > value is not in agreement with three body adiabatic 
calculations [8,38]. 

(ii) For the iV-body system, Eqs. (11-16), we have found a general expression for calcu- 
lating the nuclear fusion rate at thermal equilibrium, Eq.(25). 

(hi) Our results do not support the idea of averaging the fusion reaction cross-section over 
the momentum distribution, postulated in [6,7] and used in many publications [6,7,9-18]. 

I thank N.J. Giordano for providing the opportunity to finish this work. I am also grateful 
to N.J. Fisch for valuable discussions and for making me aware of Ref. [9]. 



6 



References 



1. J.C. Kimball, J. Phys. A: Math. Gen. 8, 1513 (1975). 

2. E. Daniel and S.H. Vosko, Phys. Rev. 120, 2041 (1960). 

3. V.M. Galitskii and V.V. Yakimets, Zh. Eksp. Teor. Fiz. 51, 957 (1966) [Soviet. Phys. 
JETP 24, 637 (1967)]. 

4. P. Eisenberg et al, Phys. Rev. B 6, 3671 (1972). 

5. L.P. Kadanoff and G. Baym, Quantum Statistical Mechanics ( W.A. Benjamin, New- 
York, 1962). 

6. N.L. Aleksandrov and A.N. Starostin, Zh. Eksp. Teor. Fiz. 113, 1661 (1998) [JETP 86, 
903 (1998)]. 

7. A.N. Starostin and N.L. Aleksandrov, Phys. Plasmass 5, 2127 (1998). 

8. G. Fiorentini et al., Phys. Rev. C 67, 014603 (2003). 

9. A.V. Eletskii, A.N. Starostin, and M.D. Taran, Physics - Uspekhi 48, 281 (2005). 

10. Y.E. Kim and A.L. Zubarev, Jpn. J. Appl. Phys. 45, L 552 (2006); Jpn. J. Appl. Phys. 
46, 1656 (2007). 

11. M. Goraddu et al, Eur. Phys. J. B 50, 11 (2006). 

12. M. Goraddu et al., Physica, A 340, 490 (2004). 

13. M. Goraddu et al., Physica, A 340, 496 (2004). 

14. A.N. Starostin et al., Plasma Phys. Rep. 31, 123 (2005). 

15. A.N. Starostin, V.I. Savchenko and N.I. Fisch, Phys. Lett. A 274, 64 (2000). 

16. A.N. Starostin et al., Physica A 305, 287 (2002). 

17. A.N. Starostin et al., Physica A 340, 483 (2004). 

18. V.I. Savchenko, Phys. Plasmass 8, 82 (2001). 

19. S. Engstler et al., Phys. Lett. B 202, 179 (1988). 

20. U. Greife at al., Z. Phys. A 351, 107 (1995). 

21. J. Kasagi at al., J. Phys. Soc. Jpn. 64, 3718 (1995). 

22. H. Yuki at al., J. Phys. G 23, 1459 (1997). 

23. H. Yuki et al., JETP Lett. 68,823 (1998). 

24. K. Czerski et al., Europhys. Lett. 54, 449 (2001). 

25. J. Kasagi et al., J. Phys. Soc. Jpn. 71, 2881 (2002). 

26. F. Raiola et al, Eur. Phys. J. A 13, 377 (2002). 

27. F. Raiola et al., Phys. Lett. B 547, 193 (2002). 



7 



28. J. Kasagi et al., J. Phys. Soc. Jpn. 73, 608 (2004). 

29. F. Raiola et al., Eur. Phys. J. A 19, 283 (2004). 

30. C. Rolfs, Prog. Theor. Phys. Suppl. 154, 373 (2004). 

31. F. Raiola et al., J. Phys. G 31, 1141 (2005). 

32. J. Cruz et al., Phys. Lett. B 624, 181 (2005). 

33. B.L. Altshuler et al., J. Phys. G 27, 2345 (2001). 

34. M.Yu. Kuchiev, B.L. Altshuler and V.V. Flambaum, J. Phys. G 28, 47 (2002). 

35. M.Yu. Kuchiev, B.L. Altshuler and V.V. Flambaum, Phys. Rev. C 70, 047601 (2004). 

36. B.N. Zakhariev and S.N. Sokolov, Ann. Phys. 14, 229 (1964). 

37. V.V. Flambaum and V.G. Zelevinsky, J. Phys. G 31, 355 (2005). 

38. L. Bracci et al., Nucl. Phys. A513, 316 (1990). 

39. A.L. Zubarev and M.Z. Nasirov, Sov. J. Nucl. Phys. 54, 389 (1991). 

40. A.L. Zubarev, M.Z. Nasirov and E.M. Gandyl, Sov. J. Nucl. Phys. 53, 566 (1991). 

41. J.N. Bahcall and R.M. May, Astrophysical J. 155, 501 (1969). 

42. G. Gamov and E. Teller, Phys. Rev. 53, 68 (1938). 

43. S. Ichimaru, Rev. Mod. Phys. 65, 255 (1993). 

44. J.N. Bahcall et al., Rev. Mod. Phys. 54, 767 (1982). 



8 



